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Abstract 

We study a rarefied gas, described by the Boltzmann equation, between two coaxial rotat- 
ing cylinders in the small Knudsen number regime. When the radius of the inner cylinder is 
suitably sent to infinity, the limiting evolution is expected to converge to a modified Couette 
flow which keeps memory of the vanishing curvature of the cylinders (ghost effect [20]). In the 
1-d stationary case we prove the existence of a positive isolated Z/2-solution to the Boltzmann 
equation and its convergence. This is obtained by means of a truncated bulk-boundary layer 
expansion which requires the study of a new Milne problem, and an estimate of the remainder 
based on a generalized spectral inequality. 

1 Introduction 

This is a revised version of the original paper fJQ, where a term was missing in the spectral 
inequality. It is well known (see for example [TTJ [2U1 E] an d references quoted therein) that the 
asymptotic behavior of the solutions to the Boltzmann equation, in the limit of small Knudsen 
numbers, is well approximated by the compressible Euler equation for a perfect gas, while the 
viscosity and heat conducting effects are seen as first order corrections in the Knudsen number. 
To get finite size viscosity effect is more delicate because of the von Karman relation [53] between 
the Reynolds, Knudsen and Mach numbers, Re, Kn, Ma: 

Ma ~ ReKn. 

When Kn is small, either i?e _1 or Ma have to be small. Therefore, finite viscosity effects are 
attained only if one assumes that the Mach number is of the same order as the Knudsen number. 
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With the extra assumption that density and temperature profiles differ from constants at most 
for terms of the order of the Knudsen number, it is then possible to show that the asymptotic 
behavior of the solutions to the Boltzmann equation is well approximated by the incompressible 
Navier-Stokes equations (INS), in the sense that the average velocity field, rescaled by the Mach 
number, converges to a solution u to INS. Moreover, the first order correction to the temperature 
profile converges to a solution to the heat equation with a convective term due to the rescaled 
velocity field u. Such results have been proved in several papers. An overview is provided in |17j . 
to which we refer for a partial list of references on the subject. We also stress that the asymptotic 
behavior of the solutions to the compressible Navier-Stokes-Fourier equations, in the low Mach 
number limit and with the same assumption on density and temperature, is the same. 

When the density and temperature do not satisfy the above mentioned assumptions, the Boltz- 
mann equation deviates from the compressible Navier-Stokes-Fourier equations. Such a discrep- 
ancy, called ghost effect [2Q], is the issue we want to address in this paper. 

The name is suggested by the fact that a small velocity field produces finite size modifications 
of the usual heat equation. These modifications are confirmed many by numerical experiments. 
There has been a big theoretical, numerical and experimental work on this and for the details 
we refer to [20l [21] . We remark that the time dependent equations on the torus, in the diffusive 
space-time scaling were written in [13] and analyzed in [9], where it is observed that similar 
equations had been previously considered in some special cases [TB] . Very little is known from the 
mathematical point of view, and the only rigorous result we are aware of is |10j . By using the 
techniques illustrated in the present paper, it is possible to deal with the time dependent problem 
in the torus, but that will be the subject of future works. 

In this paper instead, we want to study a different type of ghost effect, pointed out in [2"2l [211] 
as ghost effect by curvature. It consists in the following: if one looks at the Couette flow between 
two coaxial rotating cylinders in the limit when the radius of the inner cylinder goes to infinity, one 
expects to obtain, as asymptotic behavior, the standard planar Couette flow. The analysis based 
on the Boltzmann equation does not confirm this. Indeed, an extra term appears in the limiting 
equations, which is reminiscent of the original structure of the problem. Hence, an infinitesimal 
curvature produces a finite size discrepancy, the ghost effect by curvature. The extra term gives 
rise to a bifurcation of the laminar stationary solution, which is absent in the standard Couette 
flow. 

To be more specific, we look at the behavior of a rarefied gas between two coaxial cylinders of 
radius L and L + 1, described by the Boltzmann equation in the diffusive space-time scaling. In 
cylindrical coordinates (r, 9, z) € (L, L + 1) x [0, 2ir) x R it is written as 

dF dF vgdF dF v s ( OF dF\ l„ /ri ~ 

where the positive and normalized function F(r, 9, z, v r , vg, v z , f) is the probability density of par- 
ticles in cylindrical coordinates and we have denoted by e the Knudsen number. (v r , vg,v z ) are the 
components of the velocity in the local basis associated to the cylindrical coordinates and Q{f,g) 
is the Boltzmann collision integral for hard spheres: 

Q(f,9)(v) = lf dv* f dnB(n, v - v.){f'J + f'g' t - f t g - g,f}, (1.2) 

1 JR 3 JS 2 

with /', fi, f, /» standing for /(«'), /«), f(v), /(«,) respectively, S 2 = {n E R 3 \n 2 = 1}, B is the 
differential cross section 2B(n, V) = \V ■ n\ corresponding to hard spheres, and v, and v',v' t are 
post-collisional and pre-collisional velocities in an elastic collision: 

v' = v — n(v — w*) • n, = + n(v — u*) • n. (1-3) 
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The use of cylindrical coordinates produces a force-like term depending on the velocity, which 
will be referred below as centrifugal force. 

We will look at the above equation in the planar limit, where one takes the radius of the inner 
cylinder L to infinity. A convenient change of variables is the following: 

y + 7r 

r = L+'— — ; LB = -x; v r = v y ; v e = -v x , 
Ztt 

while the z variable is unchanged. For simplicity, we assume the distribution F invariant under 
rotations around the axis of the cylinders, so we drop the dependence on 9. Moreover, we scale 
1/L proportionally to the inverse of the square of the Knudsen number: 

1 e 2 

with a constant c, related to the curvature, which will be specified below. With these assumptions 
the equation becomes 

dF 8F dF e 2 . . ( OF 8F\ 1 fri . 

v y^ + v ^ + ^ y)v *{ v ^ y ~ v ^) = -e Q{F > n (L5) 

with 

27T+ + tt) 

The variable y varies between — ir, corresponding to the inner cylinder, and n corresponding to 
the outer cylinder. 

The boundary conditions on the two cylinders are assumed to be given by the diffuse reflection 
condition, meaning that the distribution of the incoming particles is Maxwellian: 

F(-n,z,v,t) = a-{F)M_, v y > 0, 

(1.7) 

F(n,z,v,t) = a+(F)M+, v y < 0. 
We use the following notation: for p > 0, T > and u € K 3 

\v - u\ 2 

is the Maxwellian with density p, temperature T and mean velocity u. In this paper we consider the 
two cylinders at the same temperature T = 1 and rotating with velocities U- and U+. Therefore 

M±=M(V2^,l,(U±,0,0);v), (1.9) 

where the density has been fixed so that the normalization condition 



/ 



dv\v y \M± = 1 (1.10) 

is satisfied. The constants a±, depending on the outgoing flow at the boundaries, are determined 
by the condition of vanishing net flow in the radial direction: 

/ dvv y F(y,z,v x ,Vy,v z ,t) = 0, for y = -n, n. (1.11) 

Jr 3 
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By using (ll.7[) and (|1.10l) . one immediately gets 



a-(F) = — dvv y F(-TT,z,v x ,Vy,v z ,t), 



a+(F) = / dvvyF(-K,z,v Xl Vy,v Zl t). 

Jv y >0 



(1.12) 



We need to introduce a low Mach number assumption. Due to the particular geometry we 
consider here, we do not need that the full velocity field is small. Indeed the tangential component 
can be of order 1, while we need the radial and axial components to be of the order of the Knudsen 
number e. We will use the notation v to denote the couple (v y , v z ). Correspondingly, V = (d y , d z ). 
The Mach number assumption is therefore: 

u := / dvvF = 0(e). (1.13) 

The tangential component of the velocity, denoted by U is 0(1) with respect to the Knudsen 
number e. However, we will also need some smallness assumption on it. Therefore, we introduce 
another parameter 5, measuring the size of U : 

U := I dvv x F = 0(6). (1.14) 

The fact that 5, although small, will be chosen much larger than e, is responsible for the emergence 
of a ghost effect. In principle S might be completely independent of e, but, for technical reasons, 
we will assume, in the estimate of the remainder, a specific relation between 6 and e. Therefore, 
from now on, we replace U with 5U with U = 0(1) both in e and d. In order to get well defined 
equations when 6 — ¥ 0, we will also assume the constant c appearing in the definition of a (|1.6[) of 
order 5: 

c = SC (1.15) 

for some other constant C also of order 1. 

As usual we will look for a solution to (|1.5|l in terms of a truncated expansion in e. The collision 
term forces the lowest order of the expansion to be a local Maxwellian. In order to fulfil the low 
Mach number assumptions (|1.13l) . (|1.14[) and the boundary conditions (|1.7p . the lowest order has 
to be of the form 



Ms = M(l + *r, 1 + Sr, (8U, 0, 0); v) = ^^^3/2 ex P 



2(1 + <St) 



(l.lf 



where 

v = v-(SU,0,0) = (v x -6U,v). (1.17) 

Note that the functions r and r, representing corrections of order S to the density and temperature, 
vanish at the boundary because we are restricting ourselves to the case when the two cylinders 
are at the same temperature. On the other hand gradients of U of order d warm up the fluid in 
the bulk and may produce variations of the temperature and density of order S 2 . 
The solution is sought for in the form 

F = M s + <f> + s% (1.18) 
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where 

N 

$ = ^e n F n (1.19) 

n=l 

for a suitable choice of N, and 7 is a remainder. In the next section we will give the procedure to 
compute the functions F n 'a, which will be based on a kind of Hilbert expansion for the bulk parts 
B n of F n and a boundary layer expansion in order to restore the boundary conditions violated 
by the bulk terms. The computation of the bulk terms requires the solution of a rather complex 
system of equations for the hydrodynamical fields U , u, r, r depending on 8. They are 

V[r + t] +<5V[t"t] = 0, 

V -u + S{d t r + V • (ur)) = 0, 

(1 + 5r) (d t U + u • w) = %At7 + V • (tfcW), 

(1 + 5r)(d t u + u-\7u)+ V3> 2 - ^pU 2 e v = TjoAu (1.20) 

+V • (risVu + 5 — [ctiVt ®Vt + <j 2 W ® W] , 

3 5 - 

-(1 + 5r)d t r + -(1 + <5r)tt ■ Vr = k At + V(k«5Vt) + 5-q\VU\ 2 . 

Here ?y, k, cti and <ii are suitable transport coefficients depending on 1+St. We have set rj = 770+%, 
k = kq + ns, with 770 and kq the values corresponding to S = and 775, the differences. e y is 
the unit vector in the direction y, 7 is defined in (I2.18|) and 7 2 is an unknown pressure related to 
the almost incompressibility condition given by the second of the equations (|1.20p . 
When 5 goes to the equations take a rather simpler form: 

V • u = 0, 

dtU + u- W = 77 At/, (1.21) 
d t u + u ■ V?i + VT2 - -^pU 2 e y = t] Au, 



V[r + r] = 0, 

rftr+-fi.Vr = /soAr. (1.22) 

The equations have to be completed with Cauchy initial data and the time independent boundary 
conditions 

u(±tt,z)=0, U(±n,z) = U±, t(±tt,z) = 0. 

The first of the equations (|1.22p is a Boussinesq condition ensuring the constancy of the pressure 
to the first order in 5, while the second is just the heat equation with a convective term. Note 
that the equations (|1.2ip are decoupled from the (11.22)) and can be solved independently. 

The equations (|1.2ip are the equations for the planar Couette system with an extra term in 
the equation for u, representing the curvature ghost effect. Their linear analysis [20l [22] shows 
the presence of a bifurcation controlled by the parameters C and U±, with a stationary laminar 
solution losing its stability and bifurcating into two stable non laminar solutions. 

This paper is devoted to the analysis of the 1-d stationary laminar solution to (| 1 .5[) . The two 
dimensional case where bifurcation arises, will be presented in a forthcoming paper. 
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In Section 2 will be given a perturbative analysis in S of the system (|2.31[) in order to control 
the difference between its solutions and those to the equations (|2.32[) . 

Due to the presence of the centrifugal force, the boundary layer expansion, as in [31 E] , has 
to include that force. This requires the solution of a Milne problem with a force, which has been 
given in [12] in the presence of a potential force. The present force is different because it is not 
potential and depends on velocities. Therefore the arguments in |12j require several modifications 
which are presented in Section 3. Finally, in Section 4 we estimate the remainder. The key 
ingredient to do this is a generalized spectral inequality for a perturbed linearized Boltzmann 
operator, already used in [3J in the context of the Benard problem. In the inequality given in [3] 
a term is missing. To take it into account we change the scaling to S = js's and prove a suitably 
modified spectral inequality, incorporating part of the asymptotic expansion. The remaining part 
of the asymptotic expansion is of order S 2 and easy to handle. The main result of this paper is 
summarized in the following 

Theorem 1.1. Assume 5 = 723 . Then, for e and 7 small enough, there exists a non-negative, 
isolated Li-solution to the "problem 

with diffuse reflection boundary conditions 

F(^tt,v) = tM^ dvv y F{TTr,v), ±v y > 0. (1.24) 

JVy<0 

Moreover, for q = 2 and q = 00, 

|| [J , -M(1,1,(517 > > 0)]M- 1 \\ q ,2<c^, (1.25) 
where U is the unique solution to \2. 32}) and || • || 9 .2 is defined in l{4-21\ ). 
Remark 1.1. The proof of the theorem shows that the rest term is of order s% in L^. 



2 Expansions 

In this section we show how to compute the contributions F n for n = 1,...,N in (|1.19[) . A 
modified Hilbert expansion is used to compute the bulk terms. Since they violate the boundary 
conditions, we introduce boundary layer corrections essentially supported in thin layers (of size of 
the order of e) near the inner and outer cylinders. Therefore, F n is written as follows: 

F n = B n + b+ + b~, (2.1) 

where B n is a smooth function of y, while 6„ are smooth exponentially fast decaying functions of 
the rescaled variables Y ± — e _1 (7r qp y), so that they are exponentially small away from ±7r. 

2.1 The bulk expansion 

In order to compute the expression of the B n , we substitute (I1.18[) in (|1.23[) . We ignore the terms 
b^ , because they are assumed exponentially small, and equate terms with the same power of e up 



G 



to the order N. We use the short notation 

£ s f =2Q(M 5 ,f); (2.2) 

Moreover, since in this subsection the parameter S is kept fixed, we omit the index S when there 
is no ambiguity. We get the following conditions: 

LBi = VydyM, (2.4) 
LB 2 =v y d y B l -Q(B l ,B l ), (2.5) 

£B 3 = v y d y B 2 + J^X(M) - 2Q(B 1 ,B 2 ), (2.6) 

^ n— 3 

£B„ = ^.i + ^j^^Vs-k) 

h=0 

Q(B h ,B k ), n = 4,...,N. (2.7) 

h,k>l,h-\-k— n 

In (12. 7p So = M and <r^ are the coefficients of the £-power series expansion of a(y): 

oo 
h=0 

The appropriate functional space to solve the above equations is the Hilbert space JC of the real 
measurable functions on the velocity space K 3 , equipped with the inner product 

[f,g) = f dvM-\v)f{v)g{v). (2.8) 

The operator £ is defined in a suitable dense submanifold T>£ of !H and satisfies the following 
properties: 

LI) £ is symmetric and non positive: (/,£<?) — (<?,£/); (/,£/) < 0. 
L2) £ has a 5-dimensional null space spanned by the collision invariants: 

Null£ =span{V'o,---,V'4}, ( 2 -9) 
with ipp — XpM, /3 = 0, . . . , 4 and 

M 2 

xo = !; xi = ^; X2 = « y ; X3 = « z ; X4 = — ■ (2.10) 

The orthogonal projector on Null £ is denoted P, while P- 1 = 1 — P denotes the projector 
on the orthogonal complement of Null £ in !H. 

L3) The range of £ is orthogonal to Null£: (ip a ,Lf) = for any a = 0,...,4 and for any 
L4) The following decomposition holds: 

£f = -vf + Xf (2.11) 
where % is a compact operator and v a smooth function such that 

i>o(i + M) < K u ) < M 1 + \ v \) for a11 v e k 3 . (2.12) 
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L5) If g G P^-H then, by the Fredholm alternative theorem and L4), there is a unique solution 
in P^-yC to the equation 

£/ = <?, (2.13) 
which, with a slight abuse of notation, we denote by £>~ 1 g: 

/ = £"V (2.14) 

Any solution in !H of (|2.13[) can be written as 

f=L- 1 g + f (2.15) 

with fe Null£. 
L6) Spectral inequality: there is a constant c > such that 

(/,£/)< -c(P x i>P x /). (2.16) 

By L3) in order to solve (|2.4j) we need to impose the compatibility condition P(v ■ VAf ) = 0. 
It is immediate to check that this is true if and only if 

d y y = 0, (2.17) 

where 

9 = pT, p=l + 5r, T=1 + St. (2.18) 

This is just the second of (|2.31l) . 
If this condition is satisfied, then 

v y dyM = 6(%d y U + <Bd y T), (2.19) 

where 



<8 = v x v y M, 



~ v 2 -5T ir 



(2.20) 



are in P !K. We define 21 and *B as the solutions in P^-'K of the equations 

£21 = at; £*B = 03. (2.21) 

Therefore, by L5) 

2l = £ _1 2l, Q3 = £~ 1$ 8. (2.22) 

Henceforth, 

B 1 = 5(Kd y U + Qld y r) +M (^ + ^ + ! ^=J^i) := Bi+B\. : (2.23) 

with px, T\ and u to be determined by the other conditions. 

In order to solve (|2.5[) . we need to impose the orthogonality of the right hand side of (|2.6[) to 
the null space of £. 



s 



A standard computation shows that this is equivalent to the conditions 




(2.24) 



with 



r] = - / cfoQ3<B, re 
Jm 3 

Then we can compute the part of B2 as before: 




diMSL 



B 2 = -L- 1 Q(B 1 ,B 1 )+L- 1 P x (v-VB 1 )+B, 



(2.25) 



with 




v ■ U2 



T 



+ 



v 2 -3T 
2T 2 



) 



(2.26) 



and P2, U2 and t 2 to be determined. 

The same procedure is applied to solve the equation for B3 , (12.61) . The compatibility condition 

is 



Here ai are some suitable transport coefficients of higher order whose explicit expression is given 
for example in |20j . 

Note that the term in U 2 in equation (|2.28j) derives from the contribution 
/ dvv y r H{M) — 0{5 2 ). The result in (|2.28p is independent of 5 because of the scaling (|1.15p and 
hence it persists in the limit 5 — > 0. 

The equations written so far are just the system (12 .3 1[) . They represent a system in the 
unknown functions U, u y , r and r, which does not include any of the extra funtions p\, T\ etc, 
which also have to satisfy some extra conditions, for example the Boussinesq condition to the first 
order in e. 

We do not give the explicit conditions which follow in a rather standard way. It is clear that 
the above procedure can be continued to any specific order. In this paper it will be truncated at 
N = 5. Note that at each step the solution is given up to the choice of five arbitrary functions 
which are fixed in the subsequent steps. In particular the last term of the truncated expansion, 
-B/v, is determined up to its hydrodynamic part which is arbitrary. We will take advantage of this 
arbitrariness when dealing with the equation for the remainder. 

2.2 The boundary layer expansion 

We need to include in our scheme a boundary layer expansion, because the bulk terms Bi do 
not satisfy the diffuse reflection boundary condition. For example, it is immediate to check from 
equation (12.23)) that Bi cannot be proportional to the Maxwellian M± at the boundaries. Therefore 
we introduce the corrective terms , with a fast dependence on y, so that they are sensibly 




(2.27) 



This implies several conditions. We write explicitly only the following: 




(2.28) 



9 



different from oniy close to the boundary. To achieve this, bf is assumed to be a smooth 
function of the variable Y^ 1 — e^ 1 (n =p y). We define bf as the solution to the following equation: 

v vg^± + cW^( £Y± - = + C±bf, (2.29) 

with £i g — 2Q(M±, g), M± — M(l, 1, (5U±, 0, 0); v). Here a = <rip, with <p a smooth cutoff 
function 

'i ye[o,c], 

y>2C 

for some C > and with uniformly bounded derivatives. The operator is defined in the 
same way as ^ , but we replace the hard spheres collision cross section B(n, V) with "&B(n, V) 2 . 
The reason for introducing this unphysical operator is due to a technical difficulty which will be 
discussed in the next section. The parameter ■& is chosen as i9 = c % gi , and the contribution 
from £f, which should not be there, will be subtracted in the next order of the boundary layer 
expansion. We prescribe vanishing mass flux at the boundary: 



dvv y bf(0,v) = 0. 

This equation has to be solved with prescribed incoming data at F ± = 0: 

bf(0,v) = hf(v), ioiv y >0. 

The incoming boundary data are chosen in such a way to compensate the fact that Bf is not 
proportional to a Maxwellian at the boundary: h ± (v) = —B^-(±tt). The solution to the Milne 
problem in general has a finite, but not vanishing limit at infinity, achieved exponentially fast. 
Let it be denoted by bf x . It belongs to the null space of L^. The non vanishing of bf ^ is not 
good to our purposes because this contributes to the solution in the bulk. Therefore we define 
b 1 =b 1 —b loo . In this way we ensure the decay at infinity, but b 1 do not satisfy any more the 

equation (|2.29j) , because a term of the form c % gi a(^(eY ± — ir))'N(bf 00 ) appears in the right hand 
side. We will compensate it with a term in the next order of the boundary layer expansion. We 
are not yet done, because the boundary value of Ff 1 would still be incorrect for the term — bf ^ 

which is not Maxwellian. However, we have not yet fixed the boundary values of b\ and we can 
use them to compensate it, since it is in the null space of . Finally, note that, on each boundary, 
the boundary value correction due to the other boundary is not zero, but exponentially small in 
£ _1 . This will be compensated in the remainder. In conclusion fx = B\ + b\ + b\~ satisfies the 
diffuse reflection boundary conditions up to terms ^f exponentially small in e _1 . 

The equation (|2.29l) is a special case of the Milne problem we discuss in the next section. We 
note however that in the standard Milne problem the second term in the left hand side of (|2.29p 
is absent. When in the equation there is a force term, as in the present case and in [THl [21 H], 
although very small, the lack of regularity in the velocity of the solution to the Milne problem 
for v y = (the derivative d Vy bf does not exist for v y — at the boundary), does not allow us to 
include them in higher order terms of the expansion. Indeed we need to keep it in (|2.29j) which 
will be solved in a suitably weak sense, because we cannot afford to have any v v derivative of bf 
present in the expansion. This problem was already present in the case of the Benard problem 
where the force derives from a potential and the solution is given in |12j . 

The corrections to B n , for n > 1 will be b^ solving a similar equation: 

v Vd^E+ c^^( £Y± ~ = ^bt + + St (2.30) 
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with prescribed incoming data at Y ± = 0: 

vanishing mass flux at the boundary: 

m n — / dvv y bn(0,v) = 



and then define — — b^ oc . The incoming boundary data are chosen in such a way that 
/„ = B n + b„ + b~ satisfies the diffuse reflection boundary conditions up to terms exponentially 
small in e -1 , The source term 5^ has the following form, for n > 1: 

h,k>l ,h-\-k=n 
-2 _2 



C 2 S 2 Vtv „ K n -xj ^ 

with a c = er(l — (ys) and 6q = and the property J dvS^(y, v) — 0. The existence of the solutions 
to (|2.29j) . (|2.30[) with the prescribed conditions follows from Theorem 13. 1[ given in next section, 
via a procedure which is the same used in [TH [HI [2j [3] . We do not repeat it here and refer to 
those papers for details. 

2.3 Hydrodynamical expansion 

In this section we compare the solution of the stationary 1-d equations for 5 > with the solution 
of the limiting equations for 6 = 0. The former are 

dy(pUy) = 0, 

d y [(l + 6r)(l + St)} =0, 

pUydyU = dy(T}dyU), 

5 

-pUyOyT = dy{KdyT) + S^Byll) 2 , (2-31) 
pUydyUy + 8y'J > 2 - ~Q% = ^0 9y Uy 



+d y (vsd y uy + 5 -[a 1 {dyT) 2 +a 2 {dyU) 2 }^j , 



with p = 1 + Sr and the boundary conditions 

u y (±TT) = 0; U x {±tt) = U ± , t(±tt) = r(±7r) = 0. 
In the limit of vanishing 6 the velocity u y and r are identically zero and 



riodiU = 0, (2.32) 



v 

1 



dyV2 - ~^PU 2 = 

whose solution is the laminar field U = U- + (3{y + 7r), with /3 = (2ir)~ 1 (U + — £/_). 
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The first equation in (|2.3ip implies, by using the boundary conditions for u y , pu y — 0. Since 
p = 1 + Sr for 5 small is strictly larger than zero, it has to be u y — 0. The equations reduce to 

d y ( v d y u) = o, 

dy(Kd v T) + 6(r]d y U) 2 = 0, (2.33) 
d y y2-^ P U 2 = ^d y (a 1 (d y T) 2 ). 

The first two equations decouple from the third and can be solved to find U and r. Then the last 
one gives J^. 

We define U = U-U. We have 

dyirjdyU) + [SdyT] = 0, 

dy(Kd v T) + 6 v {f3 + dyll) 2 = 0, (2.34) 
U(±n) = 0; t(±tt) = 0. 

The functions r\ and k are smooth functions of the temperature. The solutions are constructed 
by an iterative procedure. We therefore assume that ||t"||oo < 1 • In consequence H^Hoo and ||fc||oo 
are uniformly bounded for 5 < 1. Moreover, for S sufficiently small, |] /7II 00 > '-f and Halloo > f 
with t^o and kq the values of r\ and K for r = 0. Finally, note that d y r\ = SrfdyT with 77' = ^ 
uniformly bounded. By multiplying the first of (|2.34p by U and the second by r and integrating 
in y, after an integration by parts we get: 

f ||^L>|| 2 <c5||L>||||a y r||< C <5||aC/||||9,r||, 

yllVII 2 < ^IkHooCS 2 + ||a(7|| 2 ) < C <5||ar||(/3 2 + ||9(7|| 2 ). 

Therefore 

||^||<c5||^r||, 

||^r|| < cS(p + \\dyUf). 

The above inequalities imply that ||<?j,E/|| and ||9j,r|j are 0(5). In particular H^rHoo < 1 for 5 
sufficiently small. We omit the proof of the convergence of the approximating sequence, which 
follows along the same lines. 

We conclude with the following 

Theorem 2.1. If 5 is sufficiently small, the equations i2.31\) have a unique C°°(— tt,tt) stationary 
solution which differs from the laminar solution (namely U — U(y), r = 0) for 0(6). 

Call (p s , U s , u y , t s ) the solution of ([2~BTt (remember u 5 y = 0) and (1, U, 0, 1) the solution of the 
equations for 5 = with U the unique solution of (|2.32|) . 
Let Ms = M(p s , (SU 5 , 0), t s ) and M$ = M(l, (SU, 0), 1). Then 

Corollary 2.2. For q = 2, 00 

1 1 M~ 1 [Ms - M 5 °] \\ q ,2<c5 2 , 1 1 M" 1 [M° - M] \\ q , 2 < c5. 
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2.4 Estimates for the expansion 

The properties of F n 's constructed in this section are summarized in the following theorem: 

Theorem 2.3. The functions F n , n — 1,...,5 and ^ n ,e can be determined so as to satisfy the 
boundary conditions 



F n (=F7r,w) = M T (v) / \wy\[F n (Tn,w) - ^ n ^(T7r,w)]dw 
+# n , e (=F7r,u), v z ^ 0, 

and the normalization condition / R3x j_ T ^ dvdyF n = 0, so that the asymptotic expansion in e for 
the stationary problem \1.23)) , truncated to the order 5 is given by 

5 

$ = 5>"F„(y,«). 

The functions F n 's satisfy the conditions (here Q = (1 + 

|| QM- l F n || 2 , 2 < oo, || QM- l F n || 00)2 < oo , n = l,...,5, 

for any j. The functions ^S> n ,e ar e such that \\^n,e\\q,2,^, <7 = 2,oo, are exponentially small as 
e — > and J R3 dvvy^n.e = 0. 

Proof. It follows as in [21 [3] , using the Theorem 13.11 in next section. □ 

3 Milne problem 

In this section we deal with the Milne problem 

v y ^+Gu(Y)Ng = Lg + S, 

g(0,v) = h(v), v y >0, (3.1) 
dvv y M(v)g(0,v) = 0, 

dvS(y, v)M = 

for S and h prescribed. Here M is the Maxwellian with T = 1, p = 1 and u = (il, 0,0); in this 
section we adopt the notation 

Lf = 2M- 1 [L(Mf) + L#(Mf)], Nf = M _1 W(M/). 

and oj(Y) a compactly supported smooth function. Moreover, we assume that there is c > such 
that 

/ dvv y M(v)h 2 (v) < c. (3.2) 

JVy>0 

Note the explicit expression of Nf: 

df Of 

N f = v l^ v x v y - iiv x Vyf. (3.3) 

ov y dv x 
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The results will be applied with G = S 2 C2 > # > Gil, lj = a and it = <5f/±. The procedure 
is the same used in |12) : we construct the solution in a slab of size £ with reflecting boundary 
condition g(£,Rv) = g(£,v), Rv — {v Xl —v y ,v z ) and obtain estimates uniform in £, then we take 
the limit I — > oo. As in [12], the main point is to discuss the case S — 0. The assumptions on 
S will be given in Theorem 13.11 below. We only point out the differences with [12] . We write 
g = q + w with q £ Null L and (q 7 w) = 0. We set xo — 1, Xi = v x — U, X2 = %, X3 = 
X4 = — il) 2 + «y + v*], and q = J2a=o b a (Y)Xa- Moreover, ( • , • ) denotes the inner product 

on L 2 {R 3 ,Mdv). 

Note that 62 = 0. Indeed, by multiplying (|3.1[) by M and integrating in dv, we get 



_d_ 

dY 



(v y ,g) + Gu{v y ,g) = 0, 



because 

(l,Ng) = / d«sMt> v . 

Therefore, with 0(F) = dY'uj(Y'), we have 

W = (v y ,g)(Y) = exp{-G[fi(£) - Q(Y)]}(v y ,g)(t) = 

because (v y ,g)(£) vanishes by the reflecting boundary conditions at Y = I. 

As a consequence, q — X^a^2 b a (Y)xa- Moreover, (VyXm Xp) — f° r a iP 7^ 2. Therefore 
(v y q,q) = 0. 

Set I a = (v y x ai g) = (v y Xa,w) for a^2. The functions / Q satisfy the following equations: 
with 

§a,/9 = $a,!3 + (<^a,l — iWa,4)<5l,/3, 

for a, (3 ^ 2. Indeed, 

{Xa,Nxp) = j ' dvMx a (ilv x v y (Sp A - 1) - v x v v Sp t i ) , 
which is odd in Wj, for a^2. On the other hand, 

(x a ,Nw) = dv v y MxaW + Mv y v x w(S a> i -il<5 Qi4 ) 



Therefore, with / = {I a )a^=2, we have 

i(y) = j(^exp{(n(o-n(y))s}, 

which implies 

I{Y) = 

by the reflection boundary condition at Y = I. 

We take the inner product of the first equation in p. II) with g. By using 

(g,N(g)) = ^(v y g,g)- ^ii(v x v y g,g), 
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and so we obtain: 



^-j^(v y g,g) + ^Guj(v y g,g) = {w,Lw) + ^GUu(v x v y g, g) 

— (w, Lw) + —Gikj(v x v y w, w) + Giku(v x v y w, q), (3.4) 

because (v x v y q, q) = 0. The last but one term in the above Green identity is handled by adding in 
the cross section in the linearized Boltzmann operator the unphysical term dB( • , • ) 2 . Indeed, in 
this case the inequality (|2.16j) holds with v replaced by v + dv, with v satisfying the inequalities 
(|2,12|1 and v such that 

£ (i + M) 2 <v{v)< £i(i + M) 2 , Del 3 (3.5) 

for suitable constants z>o and V\. This allows to control the term ^y^(v x v v w, w) as long as we 
have # > 

So we have only to worry about the term uGU(v x v y q, w) for which we use the bound 

uGtt\(v x v y q,w)\ < IujGH\\w\\ 2 + ^GU\\q\\ 2 . 



We set A = (v y g,g). Note that A(0) < c by (|3.2j) . We need upper and lower bounds on 71(0). 
We can write 



G 



G, 



A{0) =A(£)exp{-(n(e)}+ I dY exp{-(Cl(£) - Q(Y)} 

- (w, (L + ~UGu;v x Vy)w)(Y) - Giko(Y)(v x v y q, w)] . (3.6) 

By the reflecting boundary conditions, A(£) — 0. Moreover, (w, (L + ^Gikuv x v y )w) > —(w,i/w) 
for $ > |C?It. So we only have to estimate the last term. We use the following bound proved later: 



||q(y)|| < y/ C +\A{0)\+c\\y/Uw\\(Y)+ / dY'\\yfiw\\{Y' 



(3.7) 



We use the bound 









/ ujdY 


f Y \\w\\(Y')dY' 


< [ uj(Y)YdY 1 


'o 


Jo 


Jo Jo 



<c \\w\\ 2 {Y')dY'. (3. 



Then we plug (|3.7[) in p.6p and use the spectral inequality, to get the following bound for |.A(0)|, 
using the fact that uj is compactly supported, 

\A(0)\<G(c + c\A(0)\), 

which implies the bound on |.A(0)| for G sufficiently small. Using this one can conclude that 



dY\\w(Y)\\ 2 < c 



(3.9) 



uniformly in £. 
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We need to prove the bound (|3.7[) . Let (3 a — (v y x a q) for a^2. Since f3 a = X) 7 ^2 -A-a^b^, 
with A a-1 = (VyXa,X-y) a positive non singular matrix, to bound /3 a is equivalent to estimate \\q\\. 
The equation for j3 a is obtained by taking the inner product of the first equation in (|3.1[) with 
VyXa- The result is 

df3 a 

— y = Gu) y] l B a , 1 /3 1 + T) a + (vyX a ,Lw), 



with 



and 



d 

■gy(VyXc>,w) - Gu(v y X a , Nw) 



Geo 



S an (l +il(5 a ,i - US aA ) - (vlxa,Xi) + (<W - ^ a ^){v 2 x ,vl)A a - 1 8 



7 0l,7 



An integration by parts shows that \(v y Xai Nw)\ < c||w||. The rest of the argument is as in |12j . 
The only difference is in the estimate of p a (0). We have 



\(v 2 y x a ,g)(0)\<(\v y \g,g) 1/2 (\v y \ 3 ,\ Xa \). 



(\v y \g,g)(Q) = 



Vyh — 



V y >0 



v y g = 2 



v v <0 



Vyh 



A(0). 



v v >0 



By using (JS^J) we then get ([STTjl . 

To get estimates uniform in I also for q we take the scalar product of the first equation in 
(|3.1|) and L~ l (x a v y ). The term on the right hand side is (L~ l {x a v y ),Lw), which is zero by the 
orthogonality property. We get then an equation for 9 Q = {v y L~ l (x a v y ), q) whose solution is 

6(F) = e-fo^( a K9a-^ L -i^ )iVygm {L -i {Vy ^ VyW){Y ) 

+ f dte-^ dsGu ^ 9Q ~^D(t) 
Jo 

where 9 and Q are suitable matrix, with Q = Qo, 7 6 7 , Q invertible, and 



D a (Y) = -Gu(Y) {L- l {v v x a ),Nw) 



d_ 



{L 1 (v y Xa),V y w) 



By the Schwartz inequality, boundedness of L" 1 and the fact that ui has compact support, the 
last integral is finite. Moreover, || w || vanishes at infinity and the first term on the right hand 
side is finite. This implies that there exists a finite limit at infinity, 0°°, of 8 and there is Yq such 
that for Y > Yq we have 

|e(r)-e°°| 2 < c || w || 2 (y). 

By the argument in [12] we get also for Y > Yq 

\b{Y)-b°°\ 2 <c\\w\\ 2 {Y). 

The other arguments in [12) can be adapted in a similar way. Using above estimates, the exponen- 
tial decay of ||tw|| and \b a — b™\ is established. The properties of the derivatives are also obtained 
with the method presented there, with a minor modification due to the special structure of the 
force in this case. Indeed, the v x and v y derivatives of g satisfy in this case a coupled system of 
equations. Therefore, they are both controlled only away from the boundary. 
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To state the final theorem, we define the norms 



dvM(v) / aT|/(l»| 9 



(3.10) 



for 9 > 0. 
Theorem 3.1. 

1 ) Suppose that for some (3 > 

II e? Y S || 2 , 2 ,o< oo, || e^S || OOj3 , < oo . 

J7ien £/iere is a unique solution g e L 2 (1R 3 , LooQR + )) H £2(K + x K 3 ) to t/ie MiZne problem \3.1\) . 
Moreover there exist constants c and d such that g verifies the conditions: 



Mgdv = 0, goo £ Null L, 
|| e^ Y (g(Y,v)- goo (v)) || 2 , 2 < c, || e^' y ( ff (Y, «) - ffoo («)) || 00j2 < 

/or any f3' < d . 

2) Suppose that for £ > 1, 9 > and i = 1, . . . , 3 



/ Af 


d e h 


< oo, 




< oo 


IVy>0 


dvj 




dvf 


2,2,9 



< oo. 



oo.2,e 



T/ien i/iere are constants d and q swc/i £/ia£ 



2,2,e 



oo,2,e 



for any /?' < d . 

4 The remainder 

4.1 Equation for the remainder 

It is immediate to check that the remainder 3i in (j!.18[) has to satisfy the following equation: 



d3l 



dJl 851 



l r 



dy S 2 C 
with the inhomogeneous term A given by 

A = s h+k - 2 Q(F h , F k ) - e N -\d y B N - 

h,k>lM+k>N 

2 oo N 2 N-3 

-<k* E E^ h) ^)-^5 

/s=iV-2n=0 h=0 



£,531 + 2Q(*,31) + 0(31,31) + A, (4.1) 



JV 



=W-2 



E 



, = 



C 2 S 2 



(4.2) 
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In this section we need to assume S = 763 with 7 sufficiently small. 

We recall that Ms = M (1 + 5r,l + St, (SU, 0, 0); v), where (r, r, (U, 0, 0)) is the solution of ([23Tj) . 
We keep now the dependence on 5 and denote by M the standard Maxwellian. We write 31 — MR 
and denote by L the operator 

Lf = M- 1 L Mf. (4.3) 

We use the Hilbert space H of the measurable functions on the velocity space K 3 with inner 
product 

(/,«?) = / dvM(v)f(v)g(v). (4.4) 



Note that this involves M, while the one considered in Section 2, (|2.8I) involved M7 . The operator 
L has the same properties already mentioned for L$ and we do not repeat them here. We just 
note that the null space of L is given by the span of the x/s defined in (|2.10p . We still denote the 
orthogonal projector on this kernel by P, and on its orthogonal complement by P = 1 — P. 
Moreover, we write 

M^LsMf = Lf + M _1 [£a - &o]Mf. (4.5) 



We also use the notation 
Note that 



J(f,g) = M~ 1 Q(Mf,Mg). (4.6) 



M s -M 1 + Sr , 1 „ s , 9l u 2 

-V- = (TT^ exp( -^TT^) ^ - ^ + ^ + 2 } - 1 

l$Tv 2 + SUv x + Sr - ^St + 0(S 2 ). 



By (2.1) and (2.23) 



We define 



M M 



Ms-M_ e_F\ = 1 2 3 Si 2 4 

M M 2 2 M 

Hence we have 

VK = T^i + 0( 7 M), 



where 



1 3 t> 2 — 3 

Wi = S(-tv 2 + Uv x + r - -t) + e(pi + —jj—n) 



In the present one dimensional context the u-part of B\ is zero, so 



2J(W,Pf) = 2J(W 1 ,Pf) + 2J(W-W 1 ,Pf) 

= -L(W 1 Pf) + 2J(W-W 1 ,Pf) 

= L[{-5Uv x - 1 -5tv 2 -£ v 2 ri )P/] + 0(7 2 e^) (4.7) 

:= £( a P/) + 0( 7 M), 

1 £ 

a = —6Uv x ~ —Stv — —v n 



18 



Here the property 2J(Pf,Pg) = —L(PfPg) has been used. The following operator will play 
a major role: 

Ljf :=Lf + L(aPf). (4.8) 

But a = 0(6 +e), and hence the operator L(a-) is 0(<5+e). Denote by Pj the orthogonal projection 
on Kern Lj, where Kern Lj = spanjxj — L~ 1 L(aXj), < j < 4} as in [3], Section 2. It holds that 
Pj = P- (I- P)aP. 

Theorem 4.1 (Spectral gap property of Lj). There is a constant c, such that for any function f 
in H , 

-((I + Pa)Ljf, f) > c(v(I - P)((I + aP)f, (I - P)(I + aP)f). (4.9) 
Proof . First, (I + Pa)Lj is self-adjoint for the scalar product (/, g) := J M f(v)g(v)dv. Indeed, 

M((I + Pa)Ljf(v)g(v)dv = J ML((I + aP)f)(v)(I + aP)g(v)dv, 

and L is self- adjoint for (.,.). The result follows from the spectral inequality for L. □ 

Constants which, independently of the parameter e, can be made sufficiently small for the pur- 
poses of the proofs, will generically be denoted by r\. 

Using these notation, equation (14.11) for the remainder becomes 

OR e 2 . , / dR OR 

r / \ ~\ 

i 



5 F< 



LR + 2 J(W, R) + 2 J ^2e j -^,R 



M 

2 



J(R, R)+A, 



or 



dR e 2 , . ( dR 8R\ L,R H X R T/rt „, „ . . 

+ x27^ ff (»K W^r-^^r )=^- + — + J(R,R)+A. (4.io) 



dy 6 2 C 2 \ dv y dv x 

Here A = M~ X A and 



Hi/ = 2J(W, (1 - P)f) + 2J(^£^, /) + 2J(W - W u Pf). (4.11) 



2 M 



Note that, thanks to the arbitrariness left in the determination of the Fjv, we can assume that 

PA = 0. (4.12) 
The equation for the remainder has to be solved with the boundary conditions: 

R(-n,v) = a-(R)M~ 1 M- - -#(-7r,«), v y > 

(4.13) 

R(n,v) = a + (R)M~ 1 M + - -*(tt,w), v y < 0, 
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where 

n 

M-f(TTT,z,v,t) = ^2e n V n , E (Tir,z,v,t) (4.14) 

71=1 

which are exponentially small in e -1 because of Theorem 12.31 and 

a-(R) = - I dvv y M[R(— k,v) + -*(-■ k,v)], 



(4.15) 



a+(R)= I dvv y M[R(n,v) + -V{tt,v)] 



4.2 Estimates for the remainder. 

We will proceed with the construction of the solution by iteration, based on estimates for a 
linearized problem where the non linear term J(R,R) is computed at the previous step of the 
iteration. The generic term of the iteration will then satisfy a linear equation of the type 

9R £ 2 , n , 9R d R . l rT „ TT . , „ „. 

vy— + ^7^a{y)v x (v x - v y — ) = - [LjR + H 1 (R)]+ g. (4.16) 



e 



At the n-th step of the iteration, the term g will be replaced by A + eJ(R n , R n ). Therefore 
we assume 

Pg = 0. (4.17) 

The boundary conditions are 

R(Ttt, v) = —± I {R{Tir,w) + -^{Tir,w))\wy\Mdw--y{Tir,w)), ±v y > 0. (4.18) 

M Jtw v >0 6 £ 

In order to simplify the argument we also assume that the velocity of the inner cylinder vanishes, 
so that M- is the standard Maxwellian, up to the normalization. The argument can be easily 
adapted to the general case. 

Remark 4.1. Note that the assumption (|4. 1 T[) and the diffuse reflection boundary conditions 
imply that any solution to (|4.16|) is such that 



R Vy := / dvvyRM = 0. (4.19) 
Indeed, by integrating (|4.16[) on velocities one gets: 

SyRvy + ^2°(y) R v» = 0- (4.20) 

By the boundary conditions then R Vy = 0. 

We first consider the linear problem with g given and will get L 2 bounds. This is done in two 
steps: first, we get a control of the non-hydrodynamic part in terms of the hydrodynamic one. 
Second, we get an estimate of the hydrodynamic part in terms of the non-hydrodynamic one, and 
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then a bound for R in terms of the known term g. 
The norms used below are defined as follows: 



II /Ik- = I l^dvM(v) yj^^dy\f(y,v)\ q J j . (4.21) 
Defining the ingoing velocity spaces M.± at y = =F7r as the sets v = (v Xl v v ,v z ) such that v z ^ 



||/|U,2,~=sup / dv\v y \M(v)(\f(T7r,v)\^ . (4.22) 
± \Jmi J 

We use the traces 

Note that the norm || • 1 1 2,2, is defined only for incoming velocities. In the sequel, with an abuse 

of notation we will denote by || 7"/ ||2,2,~ the || • ||2,2,~-norm of Sj~ f, where 5 is the reflection 
of the y component of the velocity. 

We now establish the lemmas which allow to bound the solution to (|4.16[) in these norms. 
• Step 1 

Recall the notation Q(v) = (1 + M) 3 , j G N. 

Lemma 4.2. Assume 5 = 7£§ , and that || (sg ||| 2 < 00 . Then, for 7 small enough and for any 
i] > 0, the solution of J^. i6'[ ) satisfies 

- || v*(I-P)(l + aP)R ||i, 2 <cU|| PR ||1 !2 (4.24) 
e L 

+e* || C| (/ - P)g |||,2 +^ II (1 + l« I) 2 * lli, 2 ,~ +^ II C|(/ - P)s lll, 2 

Proof. Define 

k(y) = exp (£ -±-^o-(q)dc^ . (4.25) 
Multiply (|4.16j) by 2i?M(l + a)fc(y) and integrate over v. We get 



d y (v y R, k(l + a)R)) = 2fc[e- 1 ((l + a)R, LjR) + e -1 ((l + a)i?, #i(i?)) 
S 

^C 2 ' 



+ {l + a)R,g)] + ^—rka I dvR 2 (5 + e)Uv x v y M + k I dvMv y (d y a)R 2 



or, with R = R^/kJyj, g = g^JHy) and * = ^./kljj), 

d y (v y (l + a)R,R) = 2[e~ 1 ((l +a)R,LjR) + £~ 1 ( J R(1 +a),H 1 (R) 



+(&(! + a), g)] j dvR 2 {S + e)Uv x v y M + / dvMv y (d y a)R 2 
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We note that, by Remark 14.11 the last term vanishes, when R is replaced by PR. Indeed, recall 
that a = —(5Uv x + \5tv 2 + ev x ll + |v 2 ti). Therefore, v y d y a is odd in v y . Since R Vy vanishes by 
Remark |4. II we conclude that 

' dvMv y {d y a){PRf = 0. 



Next, we have to look at the last term when R 2 is replaced by 2PR(1 — P)R. This is bounded by 

CS\\PR^\\v x l\\ - P)R\\ 2 2 < C(5 2 \\PR\\ 2 2 + \\v l/2 {l - P)R\\l). 

The remaining part of the last term when R 2 is replaced by ((1 — P)R) 2 can be estimated as H\ 
below. The preceding force term is of order ^ , and can be estimated in the same way. 
We next discuss the bounds for ((1 + a)R, LjR) + ((1 + a)R, H 1 (R)). First 

((1 + aP)R, LjR) < -c || i/i (I - P)(l + aP)R \\ 2 2 , 

\(a(l - P)R, LjR)\ < C( V || - P)(l + aP)R g +^||^(1 + \v\ 2 )(I - P)R\\ 2 ). 

|((1 + afaHti&ftl < C(6 3 || PR \\ 2 2 +6 \\ - P)R \\ 2 2 + 

S 2 J dvM\vf\(l- P)R\ 2 ). 

Here in e.g. the last term, after integrating with respect to y, the order of \v\ can be reduced by 
writing v = —L + K, and using a standard estimate for K (see |19j). 

(( 3+i _Kf) 2 Mdv <C J (QffMdv. 

Multiply (|4~T6| by 2RMC,jk(y) and integrate over v. We get 

^{vyRMjR)) = 2k[e- 1 {QR,LR) + e- l ^ J R,2J{W 1 R)) + (QR,g)]. 

It follows that 

e J k\v y \( J j-(R) 2 Mdv + J kvC J R 2 Mdvdy = J kQRK(R)Mdvdy 
+e J k\v y \Qj + {R) 2 Mdv + 2 J kQRJ(W, R)Mdvdy + e J kQgRMdvdy. 



Hence 



J Q +l R 2 Mdvdy <c({e + 5) J ( j+ iR 2 Mdvdy + J C j _±R 2 Mdvdy (4.26) 
+e [ k\v y \Qj+(R) 2 Mdv + e [ Q^Mdvdy) . 



In the present case j = 4, and we use (I4.26[) three times and the assumption || ||2,2 is finite. 
Moreover, 

(R(l + a),~g) = ((I P)((l + a)R), (I - P)~g) < eg || u*(I P)((l + aP)R) \\ 2 
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+ 1 e\\( i (I-P)R\\l+e3 || C l (I-P)g\\l 
Putting all the estimates together, 

+ ^ || yi(I-P)(l + aP)R\\l 2 <c[ei || C } (J - P)g |||, 2 
y |^|C47 + ( J R) 2 Md W + e5 || C§(/--P)§lla, 2 +^11^111,2]. ( 4 - 2 ?) 
where 23 := + a), R 2 {-ir, v)) - (v y (l + a), P 2 (tt, v)). 

Following the argument in [16] . one can show that the first term in 23 is bounded by 
c(2 || i/4 (J - P)(l + aP)P ||i a +t ?£ || Pi? HI 2 + £ - 3 || * HI ). 

For the second one, we shall first simplify by changing the boundary conditions at ir from diffusive 
to given ingoing data. The following lemma states the equivalence between the two problems. 

Lemma 4.3. Equation J^._?6'[ ) with the new boundary conditions 

M_ f 1 

R(-TT,v) = ^TT / (P(-7T,U;) + -*(-7T,U)))|w7 a |Mdw7 



^ J-w v >0 £ 

--V(-7r,w), v y >0, (4.28) 

R(tt,v) = —± f -^(n,w)\w y \Mdw - -V(n,v), % < 0, 
M J Wy >o £ £ 

has the same solution as problem \4.1b^ , j4-18\ l 

Proof. We start by noticing that existence and uniqueness for the new problem are classical. The 
main point is that the new boundary condition implies 

R(ir,v)v y Mdv = 0. 

Vy>0 

In fact, since J gMdv — 0, by multiplication of (|4.16l) by k(y)M and integration over v £ M. we 

get 

(k(y) J v y MRdv\ = 0, (4.29) 

which implies 

k(w) / v y MR(ir,v)dv = fc(-vr) / v y MR(-ic,v)dv. 



Hence, since J R3 v y MR(—ir, v)dv — 0, we have also J v y MR{-K, v)dv = 0. We write the integral in 
the left hand side by using the boundary condition (I4.28[) 

= / v y MR(n,v)dv = f v y M + dv f -^{-K,w))w y Mdw (4.30) 

JM 3 JVy<0 Jwy>0 £ 

VyMR(n,v)dv — f v y — '<&(ir,v) = f V y MR(w,v)dv, 

Vy>0 JVy<a S JVy>a 
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because J v <Q v y M + dv = — 1 and J R3 Vy^iw, v)dv = 0. 

Hence, J v R(TT,v)v y Mdv — 0, and this implies that the unique solution R of the new problem 
also satisfies the old boundary conditions (|4.18[) . □ 

Remark 4.2. Note that the condition J R3 R{ir, v)v y Mdv = is crucial to make this argument 
work. However, while the old boundary conditions are constructed in such a way that this is true, 
the new boundary conditions do not ensure that it is automatically satisfied. It is the conservation 
law (|4.29p which ensures the vanishing of the outgoing flux also with the new boundary conditions. 

We now return to the estimate of the second term in 23. Writing 
(1 + o)(7r) = 1 - SU+v x = (1 - \5U+v x ) 2 - \5 2 U\vl 

it follows that the outgoing part of — (v y a(n), R 2 (tt, v)) is bounded from above by 0(<5 2 ). The 
ingoing part is bounded from above by 

£ jMdv(l + 6U + \v x \)\v y \&(n,v). 

Replacing in (|4.27|) we get 



- || ^ (/ _ P)(l + aP)R ||2 < c s s || C|(/ _ P )g f 
+r/e || PR || j 2 +S 2 Ul [ Mv 2 x v y R 2 {'K 1 v)dv + [ \v y \( 4 -f + {Rf M dv 



(4.31) 



Vy>0 

-d || £s(I-P)g ||l )a +^ || * ||i, 2 ^ ]+!/(! + k, )|r,, *-(,t. r)M,lr. 



Define 



ki(y) = cxp ^3 J 



S 2 C 2 



a{q)dq 



(4.32) 



and set v x := maa;{|« x |, 1}. To estimate the outgoing integral of R in the right hand side, multiply 
(IT. 16)) by MRkiv 2 , integrate in velocity over the region v y > q, then over space using a smooth 
cut-off function \(y) which is zero close to y = —ir, and equal to one close to it, and finally over q 
for qo < q < for |go| small enough. Since we do not integrate over all v € R 3 , we cannot use the 
spectral inequality to control the terms involving Lj, but will use only the boundedness of Lj. 
Notice that norms involving \f~kR or y/kiR are equivalent to the ones for R since k(y), k\(y) are 
functions uniformly bounded in s. The notation R will be used for both. We obtain 



dq 



<lu 



v y >q 



dvMv 2 x v y R 2 {-K,v) < c\q \ - || PR \\j 



5 2 e 



(I-P)(l + aP)R\\ 2 2 . 



(4.33) 



+ - II vHl + v 2 )(I-P)R\\ 2 2:2 



° f ~ 
dq / dvdyMx'{y)vlv y R 2 (y,v). 

q Jv y >q 



The term on the l.h.s. equals 



M II «x7~-KW II 2,2— + I dq I dvv 2 x v y MR 2 (ir,v), 



'la 



q<v y <0 
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where 

-)~R(it) = R(tt, v), v y > 0. 

We move the second term in the previous expression to the r.h.s. of (|4.33[) and bound it, by 
replacing R by the ingoing boundary data, as 



dq 



0<v y <q 



M 



:VZV„dv 



Mi 



Wy >0 

2 



dww y M(R(TT,w) + -9(w,w)) 



-M~tf(7T,tf) 



<c(q )[\\ 7 -R(7r) ||L^+3 || (1 + H)*|| 



with c(qo) = o(\qo\). We replace this estimate in ()4.33p and divide both sides by \qo\, 



(4.34) 



II b*7~-RM |||^< c\ - || Pi? || 



2.2 



|| (7 - P)(l + aP)P ||| 2 +- || v-Hl P)~g || 
-~ || vHl+v 2 )(I-P)R\\i 2 }+±\\ (1 + H)§ || 



(4.35) 



It is easy to bound also 7 tt), and by similar arguments to obtain 

II 7~-R(-t) Il2,2~< c[»7 || PP ||l i2 +i || (J - P)(l + «P)P || 2 ,2 

+ - II "'Hi - P)~9 111,2 +4 II (1 + |«|)* 111,2,- ] • 



(4.36) 



This gives an estimate for the ingoing boundary term in (|4.31[) . when using the boundary condition 
in the form (l4~28)) . We also use the bound (jOS]) for || w x 7"P ||| j2 ^ in (l4~3i"j) to get, 

- e || i/i(/-P)(l + oP)fl||l ia < [£ 7 3 + ^] II ™ Hl,a 

+^ || C|(7 - P)g ||L +^ || (1 + M) 2 * ||l j2 ^ + £ l || (1 + |«|)t(J - P)g ||L . (4.37) 



□ 



• Step 2 

We will provide an a priori estimate for the hydrodynamic part of P. The equation for the 
remainder is 

dR e 2 , . , dR OR, l rr „ „ ,„ „„. 



with the boundary conditions 

M 



P(-7r,u) = I (R(-w,w) + -9(—n,w))\w y \Mdw 

-w„>Q e 



-i*(-vr,w), «j,>0, (4.39) 

R(n,v) = — i / i*(7r,iu)|w y |Mdiu- i\&(7r,v), w a < 0. 
Jw„>a £ ' £ 
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Lemma 4.4. If Pg = 0, then the solution of ^.38\) , j4-39\ ) satisfies 

II PR \\l, 2 < || !/*(/ - P)(l + aP)i? ||1 >2 + || g ||| j2 +1 || viz ||i ^ ] . (4.40) 

Proof. Let i? = 7 y R be the Fourier-transform in y of i? and ^ € Z the conjugate variable to y. 
i? satisfies the equation 



dv x 



iVytyR+^Zy \ a(y)v x (v x — - v v — ) }> = e^iLjR + HxiR)) - Vy r(-1)^ + g (4.41) 



with 



r(t;) = -R(7r, u) — i?(— 7r, w). 



(4.42) 



With the notations < R > for the 0-Fourier term, and i? := i?— < R >, we shall first give an 
estimate for R. We use the notation 



Z = s^iLjR + Hx(R))- 



6 2 C 2 ' 



3 y < a(y)v x {v x 



dR 

8V V 



OR 
' dv x 



)>+g-(-l)^v y r, 



Z = £- l (£jR + H 1 (R)) +g- (-1)*» V, 
Z' = £- 1 (lJr + HiiR)) +g,U= {ityVyY 1 . 

Let f) be the indicatrix function of the set {v € ffi. 3 ; \v y \ < a}, for some positive a to be chosen 
later. For £ y ^ 0, 

ii ptoHv,')) ii 



4 
3=0 



dvt)(v)R(€ v ,v)M(v)xj(v) 



< c ii c-M ii E ii to ii < <vs ii c- s f)P ii • 

3=0 

We use this estimate with the following choice of a, a —\\ C-sR ll^ 1 ! C-sZ' II- We also introduce 
an indicatrix function t)i with a — 5\. We fix 5\ so that cS\ << 1. Then we find from the above 
estimate that the P-part of the right-hand side, || P(f)i-R) ||, can be absorbed by || P(f)i-R) || in the 
left-hand side. The estimates hold in the same way when f)i is suitably smoothed around 
For the remaining fff^P = (1 - f))(l - f))P, we shall use that R = -UZ. Then 



ii pm c k{i v ,-)) f 

< C || Qs+2mU f\\ C-sZ' t 



pm c v y 



Sl\iy\ 2 



< 



CsZ 



I 1 1 2 



pm c v y r) 

Si\^y\ 2 



5 2 C 2 



5 2 C 2 



|6|, 



|0| 



with 



9= E / X3^ c U? y U(y)v x (v x ^ - v y ^) 

j=o ^ y x 

Mdv ( I X j^ c {R~UZ)Mdv 
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We replace a by || (~ S R || 1 || C-sZ' | in the denominator. That gives 

,2 ^,„„ d ii ii a , \\ p m c v y r) 



|2 



pp(e v , on 2 ^caic-izimc-s'ii . — ^r^- 



.2 



+*i || C_ a (J-P)fl|r) + ^a|e|. 



Hence, 



PR(Zy, • ) ll 2 < c((|| PP || + || C-.(/ - P)R IDs II C-^' II 

ii p(f)frv) ii 2 



Consequently, 



PR(Cy, • ) ll 2 < c( || C-^' II 2 + lim p l lj ) lr + II - P)P II II C-^' 



+ ||C- s (/-P)i?|| 2 j+^|6|. 
We next discuss the term s l cS |Q|. The first integral can be bounded by an integral of a product of 

2 i 

M, 1 + \£ y \, a polynomial in v, \ R | and U or U 2 . Since, by our choice of 5, we have = ^§^7, 
this integral is bounded by eic \\ R ||2- And so, 

11 ppfe, • ) n 2 < «( 11 c-sz' 11 2 + " p( ^y l|2 + 11 ( j - 11 

Therefore for £ y 7^ 0, 



*iievi 

2 



y |pi?i 2 (^, v )Md V < c (i 11 c- a («)Ljfl(e v ,o ni +^ 11 c- s («fe,') 



+ || (/ - P)P(£,, •) || 2 + " ; " + II -)H 2 J ■ (4-43) 

The error from evaluating P instead of Pj is of order S+e. We remind that the zero Fourier mode 
of PR is zero by definition. Hence, taking e small enough and summing over all 7^ £ y € Z, 
implies by Parseval that 

y (PR) 2 (y,v)Mdvdy < J z/((J - P)(l + aP)R) 2 (y,v)Mdvdy 

v' x g 2 {y, v)Mdvdy + 5 \\ PR \\l 2 ) . (4.44) 



P(f)!*) c v) 



(5! 

The r-term can be expressed from (|4.4ip at £ y = 



{OR OR 1 

^ )v ^oV y - v ^) ^ (445) 

+g(0,v) + -Ih{R)(0,v). 

e 
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Inserting this into (|4.44l) results in 

J (PR) 2 (y,v)Mdvdy < J v{{I - P){1 + aP)R) 2 {y,v)Mdvdy 

+ /^V<*.)«** + *i™ib). 



We are left with the Fourier component PR(£ y ) for £ y — 0, and have 



PR(0,v) = Y,xe(v) J dvM^ J dyR{y,v). 



We discuss each moment separately. Given two functions h(v) and /( • ,v) we use the notation 
fh{' ) '■— J dvh(v) f (■ , v). In particular, for h = Xjt 3 = 0> • • • > 4, we also use the notation fj. 

• Wj,-momcnt: 

As already noticed in Remark 4.1, multiply (j4.38[) by k(y)M and integrate over velocity. Now 
integrate over [— ir,y} and v € R 3 . Since J v y M(v)R(±ir,v)dv = 0, we have from (I4.38[) R v (0) = 
/ v y M{v)R(y 1 v)dvdy — 0. 

• w x -moment: 

We estimate the moment R Vj . v 2 (0) and then use that 

Rv x = Pv x v? — P Vn: vp 

with i?- 1 = (1 - P)R. 

Indeed, J MPRv x v 2 dv = J MRv x dv, because / Mv 2 v 2 dv = 1. Multiply equation (j4~38]) by 
Mv x v y , integrate over [y, tt] and u G R 3 and then integrate over y £ [— tt, tt]. We get 



dy / dvv x v 2 MR(y,v)dv\ < |27r / dvv x v 2 MR(rr, v)dv\ 



9 II 2,2 



(4.46) 



£ V || Pi? ||a,a +- II {I - P)(l + aP)R || 2 ,: 

£ 



In the boundary term, the integral over v y < is easy because the boundary condition in tt is 
given in terms of 'J. To control the outgoing (v y > 0) part, we multiply again equation (|4.38|) by 
2irMv x v y , integrate this time over [— 7r, tt] and v £ M. 3 ,v y > 0. We get 



1/ dvv x VyMR(n,v)dv\ < | / dvv x v?.MR(-7r, v)dv\ 



Vy>0 



(4.47) 



Vy>0 



+c 



e*l 2 || PR ha +- || (/ - P)(l + aP)R || 2 , 2 + || g || 2 , 2 



To control the second term in l.h.s. of (|4.47|) . we use the boundary condition in — tt. Finally, we 
replace the bound of the outgoing part given by equation (|4.47[) in (|4.46l) . The result is 



|iW,(o)l ; 



< c[et 7 4 || PR ||l i2 



dy J v x v 2 MR(y,v)dv 
~ || (I - P)(l + aP)R ||| 2 + || <? ||l, a +^ || (1 + M) 2 * HI, 



(4.48) 



and the same estimate holds for iRy 
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• t> 2 -moment: 

we get the same estimate for the u z -moment and the proof is analogous, the only difference being 
that we start by multiplying (|4.38p by Mv y v z . 

• Xi~ moment: 

We recall Xi — ^|/g 3 and we denote Ri(0) = f*^ XiMRdvdy. We notice that 

R v2yA (0) = ^^4(0) J vyiMdv + R^ A (0), (4.49) 

where A is the non-hydrodynamic solution to 

L(v y A) = v y (v 2 - 5). (4.50) 

To control R V 2 A (Q) we multiply (|4.38l) by Mv y A and proceed as before; first, we consider the 
integral f*^ dy J* dy' J R3 dv, then we study 2n dy J v >Q dv and take the difference. Now, the 
analogous of the second term in (I4.47[) is 



2vr / dWyAMR(-Tr,v)dv. 



Vy>0 



We use the boundary condition in — ir, and observe that J v y AMdv — by orthogonality. Since 
the integral is even in v y , also J v >Q v y AMdv = 0. 

• xo-moment: 

/2<R f 
dvv y RMdv = R + — ^=r + J dvVyR Mdv, and we already have estimated Ri(0), it is 

enough to estimate the moment R v z(0). To this end, multiply (j4.38p by v y M and integrate over 
[y, 7r] x R 3 . Since v y is in KermL, we do not get contribution from the L, Hi and g terms in the 
right hand side. Moreover, by integration by parts, there is no contribution depending on PR in 
the force term. We have 



dWyMR(y, v)dv + 2tt / dvtfMRfa v)dv + 2tt / dvvlMR{-K 1 v)dv 

Jv v <0 Jv v >0 



C^f J dvMR^-vi). 



The second term depends on the ingoing flow at 7r, which is given in terms of vp. To control the 
third term, we will as before estimate it using another equation. Multiply (|4.38[) by 2ttv v M and 
integrate over [— ir, tt] x {v £ R 3 : v y > 0}. 

dWyMR(n, v)dv\ < | / dvv 2 y MR{ — 7T, v)dv\ 

Uy>0 JVy>Q 

4 II R h,2 + II 9 h,2 +- || (I-P)(l + aP)R\\ %2 

7^ e 
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We get 

| / dvvlMR{y, v)dv\ < \2ir I v*MR(-Tr,v)dv\ 

Jv y >a 

-< I -4 || R h + || 9 h +- || (/ - P)(l + aP)R \\ 2 +- || # || 2 ~ 

7 J e e 

To estimate the first term to the right, we cannot employ the previous symmetry arguments, but 
will instead use Schwartz' inequality 



/ 



v%M\R\(-v,v)dv < 



v v >0 



1 1/2 r 



v y >0 



1 1/2 



v y MR 2 (-ir,v)dv 



v y >0 



We can now use (|4.36|) to estimate the last integral, and get 

\R vl | 2 < c[r, || PR \\l 2 +1 || (/ - P)(l + aP)R + || g || 2 2 +1 || * ||2, ^ ] . 
Collecting all the moment estimates for £ y = 0, we get 

|| PP(0) |||< c[„ || Pi? II?. 2 , +1 || (7 - P)(l + aP)R ||* a + || g ||| )2 +1 || (1 + M) 2 * ||* 2i „ ] 



and this concludes the proof of Lemma 14.41 



□ 



Combining the two steps we have proved 

Theorem 4.5. Assume 8 = £%j, Pg = 0, and \\ (l + \v\)i g) || 2 , 2 finite. Then, forj small enough, 
the solution of \4.38j) , j4.39\) satisfies 

II PR \\h< II Cffl 111,3 +|r II ||1, 2 ,^ +e-i I) C|9 ||1, 2 ) 

(4.51) 

|| !/*(/ - P)(l + aP)P ||| 2 < c(l || C2* ||I, a .~ +e f II C|3 ||1,2 +e f II Cf <? 111,2 )• 
Proof. Replace gHH) in (|4~40| . and use 7 and r\ small. □ 

The step from L 2 to is done by studying the solution along the characteristics as [2]. The 
result is 

II v^R ||^ j2 < c(l || C2P ||1,2 +e 2 || ^k 2 .g ||L,2 +^ II C2* ||I,2,~ )• 
Then, by Theorem 1431 

II ^KlR ||2o )2 < c( £ -| I) C|3 111,2 +£ 2 II C|5 ||L,2 

+ IKjff 111,2 +^IIC2* 111,2,- )■ (4-52) 

With the a priori estimates provided by Theorem 14.51 we are now in the position of proving 
the existence theorem for the remainder equation (|4.16H4.lB|) . by an iteration procedure. 
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Theorem 4.6. There exists an isolated solution in L 2 ([—ir, tt] x R 3 ; Mdvdy) to the problem J^. 10\ ), 

Proof. The remainder term i? will be obtained as the limit of the approximating sequence R n , 
where R° — and 



dR n+1 e 2 . . , <9i? n+1 <9i?" +1 . l rr n „ , , „. , n 

^-BT + 6^° iy)vAVx ^ - Vy — ] = e [LjR + HlR 



+ J(R n ,R n )+A, (4.53) 



R n+1 (-7r 7 v) (« rl,1 (-7r !U ,-) + ^( • W. . 

1 



^ -'-w v >0 : 



*(-7r,w), u a > 0, (4.54) 

i?«+ 1 ( 7r ,w) — ~tt [ -*(ir,w)\w y \Mdu>-~*(Tr,v), v y <0. 
M J Wy> Q e ' e 

Here A is of order e 4 with (1 + |u|)i.A quadratically integrable by Theoreir l2.3[ and it holds that 
R n+1 satisfies also the boundary conditions (|4.18j) . The function R 1 is solution to (|4.38M.39|) with 
g = A. Then, by using the a priori estimates of Theorem 14. 5 [ (|4.26[) . and (|4.52[) . together with 
the exponential decrease of Vt, we obtain, for some constant ci, 

II C2R 1 |U,2<c ie t, || £ 2j Ri || 2 ,2< ciet. 

By induction for e sufficiently small, 

II C 2 i?" ||oo,2< 2c l£ i, 

II C 2 (i?" +1 - R n ) h,2< c 2 e k - || ( 2 (R n - R n - r ) || 2 , 2 , n > 1, 

for some constant c 2 . Namely, 

d{R n+2 -R n+1 ) e 2 . . . d(R n+2 -R n+1 ) d(R n+2 - R n+1 ) . 
Vy dy + p&'M^ dv~ y V « dv~ x } 

= -L T (R n+2 - R n+1 ) + -H 1 (R n+2 - R n+1 ) + G n+ \ 

£ £ 
( R n+2 _ R n+l ){ _ n)V) = M± f ( jR «+2_ i? »+l)(_ 7r;W )| Wy | M _ du;j Vy >0 , 

M Jw.^O 



(R n+2 - i?" +1 )(7r, v) = 0, v y < 0. 
Here, G n+1 = (I - P)G n+1 = J(R n+1 + R n , R n+1 - R n ). It follows that 



C 2 (i?"+ 2 - R n+1 ) \\ 2 ,2< c£~i || (iG n+1 



2,2 

\ II ? { nJi+1 n)n I 

2 



< C£-i ( || ( 2 R n+1 ||oo,2 + || (2R U ||oo,2 ) II (2(R n+1 - R n ) || 2 

<c 2 £^ || ( 2 (R n+1 -R n ) || 2 , 2 . 

Consequently, 

II C 2 iT +2 h,2 <|| UR n+2 - R n+1 ) h,2 +■■■+ II UR 2 - R 1 ) Ika 

+ II C 2 i? 1 || 2 ,2< 2C1£*, 
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for e small enough. Similarly, || C2-R" +2 |U,2< 2aei . In particular (R n ) is a Cauchy sequence 
in l? M ([— 7r, 7r] 2 x M 3 ). The existence of a solution i? to (|4.10j) follows. Local uniqueness follows 
along the same path. □ 

Corollary 4.7. There exists an isolated, non-negative L^-solution F to il.2,'J\) , jl.24\ ) such that 

|| M-^F-Ms] || 2 ,2<C 7 £3. 



Proof of Theorem ] 1.1[ The non- negativity can be proved similarly to [5]-[5]. Then, by Corollary 
12.21 and 14.71 we have, for q = 2, oo 

|| M- X [f -M(1,1,(<517,0,0))] || 9 , 2 <|| M-^F-Ms] \\ q , 2 
+ || M- 1 [M S -M(1,1,{SU,0,0))] \\ q ,2<c{eS + 6 2 ) 

which implies (|1.25l) by taking into account the relation 5 = ■ye^ . □ 
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